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Abstract 

Results are provided from a viscous 
shock layer (VSL) analysis of the reentry 
flowfield around the forebody of the Japanese 
Orbital Reentry Experiment (OREX) vehicle. 
This vehicle is a 50-deg spherically blunted cone 
with a nose radius of 1.35 m and a base 
diameter of 3.4 m. Calculations are done for the 
OREX trajectory from 105 to 48.4 km altitude 
range. A 7-species chemical model is found 
adequate for the flowfield analysis. However, 
for altitudes greater than 84 km, the low- 
density effects (such as thermal non equilibrium 
and slip) are to be implemented for good 
agreement between the predictions and flight- 
inferred heat-transfer rate data. Further, at 
altitudes lower than 84 km, a finite surface 
recombination probability is to be employed in 
place of a non-catalytic surface for better 
comparison between the calculations and data. 
VSL results are also compared with the direct 
simulation Monte Carlo (DSMC) predictions at 
high altitudes (> 80 km) and the electron 
number density data for three altitudes in the 
OREX trajectory. Overall, there is a good 
comparison between the flight data and 
calculated results. With the ongoing 
refinements in data extraction procedures, the 
OREX data should prove valuable for validating 
theoretical models employed in flowfield codes 
for calculation of reacting-gas flowfields. 



Nomenclature 

d 

base diameter 

Knoo 

freestream Knudsen numer, 


= XoJd 

L 

shuttle length 

Moo 

freestream Mach number 

n 

distance normal to the body 


surface 

Pw 

surface pressure 


q w surface heating rate 

RN nose radius 

s distance along the body 

surface measured from the 
stagnation point 

T V e vibrational-electronic- 

electron excitation temperature 
T w surface temperature 

Too freestream temperature 

Voo freestream velocity 

x axial distance from stagnation 

point measured along symmetry 
axis 

Xj mole-fraction of species i 

y radial distance from symmetry 

axis 

70»7N recombination probabilities for 

atomic oxygen and nitrogen, 
respectivly 

Xoo mean free path 

T| e electron number density 

Poo freestream density 

Introduction 

For the optimum design of reusable 
space transportation vehicles, it is necessary to 
determine the reentry aerothermal environment 
of these vehicles accurately so that the weight of 
the thermal protection system may be 
minimized to increase the payload capacity. An 
accurate prediction by computational fluid 
dynamics (CFD) methods of surface heating, 
temperature, and flowfield quantities during 
reentry naturally requires accurate modeling of 
the flowfield chemistry, gas surface interaction, 
body (and shock) slip as well as thermo-chemical 
nature of the flowfield 1. Due to the extreme 
physical conditions encountered and associated 
modeling difficulties, it becomes essential to 
calibrate the CFD codes, ideally against flight 
data, for a wide range of flow conditions. The 
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Orbital Reentry Experiment (OREX) was 
conducted recently by the Japanese^ from orbit 
to create such a data base to study the reentry 
technology and establish the reliability of 
computational design tools for Earth-Space 
Round Trip (ESRT) systems. 

OREX was the first of the three flight 
experiments planned to obtain the technology 
base for the development of Japan s unmanned 
space shuttle called HOPE (H-II Orbiting 
Plane). OREX was launched into Earth Orbit by 
the H-II rocket on February 4, 1994, and was 
the first Japanese entry experiment. Two of 
many objectives of OREX were to: (i) gather 
reentry data (such as those for the aerodynamic 
and aerothermal environment and 
guidance/navigation/ control), and (ii) test the 
thermal protection system that has been 
developed for HOPE. For meeting the second 
objective, the maximum heating rate 
encountered by OREX was almost equal^ to that 
computed for HOPE. The initial reports^ 
provide an indication of the scope and quality of 
the basic data acquired. 

The inferred surface heating rate data 
for OREX are considered potentially unique for 
altitudes above about 92 km. Unlike the 
Shuttle Orbiter, which is an operational vehicle, 
OREX was enclosed in a protective fairing 
during launch thereby negates the requirement 
for waterproofing the thermal protection 
materials. Consequently, the OREX thermal 
protection system should not produce some of 
the outgasing products and the associated 
reduction^ in heating measured for the Orbiter 
during the initial phase of the heat pulse. A 
motivation, therefore, exists for the comparison 
of numerical predictions with the OREX data. 
Reference 3 presented a 7-species, one- 
temperature (1-T) analysis of the OREX 
stagnation-point heating data. The 
computations of Ref. 3 were recently redone® 
with new values for the thermal properties of 
the C/C nose cap. The low-density effects (such 
as the thermochemical nonequilibrium and slip) 
were not considered in the analyses of Refs. 3 
and 5. In the data reduction procedure for the 
electrostatic probe^, however, a two- 
temperature (2-T) calculation method® was 
employed without the surface-slip effects. 

In the present study, the OREX entry is 
analyzed by using the viscous shock-layer (VSL) 


technique of Ref. 1. This technique is used to 
calculate flowfield and surface quantities along 
the OREX forebody for altitudes of 105 km to 
48.4 km. Information concerning the surface 
temperature was obtained from Yamamoto® for 
these altitudes. Comparisons of VSL results 
with flight data include stagnation-point 
heating rates and the flowfield electron number 
densities. An assessment of the low-density 
effects (such as thermochemical nonequilibrium 
and slip) on surface heating and flowfield 
quantities is made by comparison with the 
direct simulation Monte Carlo (DSMC) results 
obtained firom the method of Bird^ and the flight 
data. 

Problem Definition and Methodology 

OREX Geometry and Sensor Locations 

The OREX vehicle is a 50-deg 
spherically blunted cone with a 1.35 m nose 
radius and a base diameter of 3.4 m, as shown 
in Fig. 1(a). It was launched into Earth Orbit 
using Japan's new H-II launch vehicle. During 
entry, aerothermal data were obtained from 
about 120 km down to about 40 km, including 
the blackout period when the maximum heating 
occurred. 

Figure 1(b) shows a schematic of the 
various sensor locations on the exposed forebody 
of OREX, which consists of thermal protection 
materials that have been developed and are 
being evaluated for use on the HOPE vehicle. 

The nose cap is a monocoque carbon-carbon 
(C/C) nose having a thickness of 4 mm. 
Surrounding the C/C hot structure are twenty- 
four C/C tiles 1.5 mm in thickness. The nose 
cap and the surrounding tiles are made from the 
same carbon-carbon mix. Ringing the C/C tiles 
are some 230 ceramic tiles made of silicon-oxide 
and aluminum-oxide fibers connected to OREX’s 
aluminum honeycomb shell. The ceramic tiles 
are 20 mm in thickness. 

A few of the locations where 
measurements were made along the OREX 
forebody are also shown in Fig. 1(b). The 
measurements for which the data are compared 
with current VSL and DSMC calculations are: 

(i) the inferred heating rates extracted® from 
the back surface (with material thickness of 
4 mm) temperature measurements made at the 
nose cap stagnation point, and (ii) the electron 
number density distribution in the boundary 


2 



layer measured by an electrostatic probe 
mounted on the conical flank before the probe 
shoulder. The electrostatic probe protrudes 
from the OREX vehicle surface 70 mm in height 
as shown in Fig. 1(c). Five semi-cylindrical 
electrodes of 0.2 mm diameter situated along 
the leading edge of the probe with 0.4 mm 
bluntness and 60° sweep angle are numbered as 
shown in this figure. Those electrodes collect 
ions at 5 vertical positions to give the ion 
density profile in the boundary layer which is 
assumed equal to the electron number density to 
be measurable between 10 1® - 10^0 
electrons/m^. The lowest altitude where the 
electron density measurement can be made 
within this range is about 75 km. The 
procedure used to reduce the probe data is 
similar to that used with the RAM-C flight 
experiment^. Further details concerning the 
OREX flight measurements can be found in 
Ref. 2, 3, and 8. 

Computational Methods 

The VSL method of Ref. 1 is used to 
analyze flowfields over the OREX forebody for 
altitudes from 105 km to 48.4 km and Mach 
numbers from 27.1 to 9.1), whereas the DSMC 
method of Bird^ is employed mostly at high 
altitudes (greater than about 80 km) in the 
present study. The VSL method can be used to 
compute^ reacting gas flows in thermochemical 
equilibrium or nonequilibrium state with and 
without body and shock-slip boundary 
conditions and with an arbitrary number of 
chemical species. For the current study, the 
number of species ranged from 5 to 11 using the 
VSL method while the DSMC method used 
herein considered only 5 species. 

Detailed description of the VSL method 
used is given in Refs. 1 and 10. Briefly, this 
method is a spatial-marching, implicit, finite- 
difference technique, which includes coupling of 
the global continuity and normal momentum 
equations. For the thermal nonequilibrium 
calculations with two- temperatures, two energy 
equations are solved for the translational- 
rotational and vibrational-electronic-electron 
temperatures. Only, the total energy equation 
is needed to be solved for the thermal 
equilibrium calculations. Further, for the 
thermal non equilibrium case, air chemistry is 
modeled by using a modified Arrhenius 
expression for the forward and backward rate 
coefficients (see Table 1 of Ref. 10). Chemical 


vibrational coupling is taken into account 
through the preferential dissociation and 
recombination model proposed by Park* \ The 
body and shock-slip boundary conditions are 
implemented at higher altitudes in the OREX 
trajectory. A description of the thermodynamic 
and transport properties, chemical kinetics 
model with chemical-vibrational coupling, 
relaxation processes for the vibrational - 
translational and electronic-electron- 
translational energies, and body and shock-slip 
boundary conditions is given in 
Refs. 1 and 10. 

The surface is assumed to be 
noncatalytic in the present study for most of the 
calculations as suggested in Ref. 5. However, a 
finite catalytic wall boundary condition ^ 
(based on the Shuttle flight heat-flux 
measurements) is also implemented, which 
impacts the heating at altitudes lower than 
about 84 km, At higher altitudes, finite wall 
catalycity does not affect heating much as 
compared to a noncatalytic surface. To obtain 
an estimate of the reduction in heating due to a 
noncatalytic (or finite catalytic) wall, some 
results are also obtained for a fully catalytic 
wall, which should be closer to the 
thermochemical equilibrium value. 

The noncontinuum 2D/axisymmetric 
DSMC method, used to compute the results 
shown here, is described in Refs. 7 and 13 and is 
briefly discussed here. In this method, the 
molecular collisions are simulated using the 
variable hard sphere (VHS) molecular model. 
The collision cross-section is a function of the 
relative energy in the collision. Parameters 
used in the present study to define the VHS 
model are a reference temperature (2880 K), 
reference diameters for each of the five species 
(3.062 x lO' 10 , 3.083 x lO* 10 , 2.297 x 10' 10 , 

2.398 x 10* 10 , and 3.065 x 10' 10 m for 02, N2, O, 
N, and NO respectively), and the temperature 
exponent (set to 0.73) in the power law for 
viscosity coefficient. Energy exchange between 
the kinetic and internal modes is controlled by 
the Larson-Borgnakke statistical model. For 
the diatomic molecules, a rotational relaxation 
collision number of 5 and a vibrational 
relaxation collision number of 50 are used. The 
gas surface interactions are modeled by 
assuming the gas molecules to reflect diffusively 
with full thermal accommodation at the 
specified surface temperature. Similar to the 
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VSL calculations, the finite catalytic wall 
boundary condition is imposed by using the 
surface recombination probabilities for atomic 
oxygen and atomic nitrogen inferred from 
Shuttle data 

Results and Discussion 

As described in the Introduction section, 
it is necessary for the computer codes with many 
models for simulating the physical processes to 
be validated against a wide range of flow 
conditions if these codes are to serve as useful 
design tools. Some of the codes show good 
accuracies against ground-test results, which, 
however, do not simulate high temperature real 
gas effects simultaneously as encountered in a 
flight. Thus, the code comparison against the 
flight data can not be over-emphasized, 
especially, for the high-energy flows. 

This section provides details about the 
flight conditions for which calculations are 
made; the sensitivity of the calculated results to 
assumptions concerning flowfield chemistry, 
wall boundary conditions and freestream 
density; flowfield and surface quantities; 
comparisons of these quantities with flight 
measured data and, at higher altitudes, with 
direct simulation Monte Carlo (DSMC) method; 
and a brief discussion of the potential 
uniqueness of the OREX heating data set. 

Flight Conditions and Numerical Parameters 

Table 1 provides the range of flight 
conditions considered, which encompass 
altitudes of 105 km to 48.40 km. Figure 2 shows 
the atmospheric (freestream) density from three 
different sources for this altitude range. The 
atmosphere values denoted as OREX are those 
given by Yamamoto®* ® for which he has 
performed Navier-Stokes calculations. These 
values are similar to those for the 1962 U.S. 
Standard Atmosphere up to an altitude of 
about 90 km. For altitudes above 90 km, the 
atmospheric properties given by Jacchia*® (for 
an exospheric temperature of 1200 K) are 
employed. Jacchia values have been used 
extensively in previous high-altitude studies 16 > 

* . These values are different than those given 
by the OREX atmosphere® and avoid a 
significant change in the slope of the density- 
versus-altitude curve around 90 km altitude (see 
Fig. 2). The density values from the Jacchia 


atmosphere are lower by about one-third from 
those given by Yamamoto (see Table 1) at higher 
altitudes. 

The stagnation-point temperatures (see 
Table 1) and surface temperature distributions 
(see Fig. 3(a)) employed in the present study are 
those obtained by Yamamoto® for altitudes of 
105 km and below using a Navier-Stokes 
solution for the flowfield coupled with a material 
response code as described in Ref. 3. The 
considerable temperature variations shown in 
Fig. 3(a) are due to variations in materials, 
material thicknesses, and material thermal 
properties. The carbon-carbon nose cap with a 
thickness of 0.4 mm extends to a wetted length 
of s = 0.942 m (measured from the stagnation 
point). Following the nose cap are the carbon- 
carbon tiles which terminate at s = 1.242 m and 
then four rings of silica tiles, having much lower 
conductivities and hence higher surface 
temperatures, extend to s = 1.982 m. The comer 
shoulder is also protected with a fifth ring of 
silica tiles. Figure 3(b) shows a comparison of 
the history of calculated values and flight data 
for the rear surface stagnation-point 
temperature. The computations of Ref. 5 for 
the surface temperature are different from those 
of Ref. 3 for two reasons: (i) thermal properties 
of the C/C nose cap are changed to new values, 
and (ii) internal emission from the C/C nose is 
assumed to be zero by the heat shield effects, as 
recommended by the industry. These two 
changes result in a C/C nose cap temperature 
history, calculated by the coupled Navier- 
Stokes/material response (CFT)-FEM) code, to 
agree better with the flight temperature history 
for the altitude range of 105 km to 45 km. 

The VSL calculations have been made 
by using a normal grid with 101 points, which 
are clustered both at the body surface and shock 
to capture large gradients in the flow properties 
there at low densities. The points are clustered 
only at the body surface at higher densities. 

The grid spacing is varied from 1x10'® m to 
lxl0~4 m to obtain a grid-independent solution 
for the altitude range given in Table 1. The 
smallest grid is used at the lowest altitude of 
48.40 km and is increased in inverse proportion 
to the density for higher altitudes. In the 
stream wise direction, a minimum grid spacing 
of 2x10*2 m is used. 

For the DSMC calculations*®, the size of 
computational cells adjacent to the body surface 
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in the direction normal to the surface is usually 
less than half of the local mean free path length, 
and for most solutions much smaller. Previous 
experience has shown that such a resolution 
a^acent to a surface provides results that are 
independent of further cell refinement. For the 
lower altitude cases, the computational cells had 
very large aspect ratios with the dimensions 
along the surface equal to many local mean free 
path lengths. 

Flowfield and Surface Quantities 

This section gives results from the 
detailed VSL calculations for flowfield and 
surface quantities at altitudes of 92.82, 88.45, 
84.01, and 59.60 km (peak-heating altitude in 
Ref. 3). These results are later compared 
extensively with the flight data for electron 
density and stagnation-point heating. 

Presented results in this section include an 
assessment of the thermal nonequilibrium (i.e. 
2-T v/s 1-T predictions) and slip effects on 
flowfield and surface quantities at these 
altitudes and on stagnation heating along the 
entire flight trajectory. A comparison of the 
VSL results is also made with the Navier-Stokes 
calculations of Yamamoto®, and the present 
direct simulation Monte Carlo computations. 
Before comparing the present results with other 
calculations and flight data, however, it is useful 
to have an indication of the sensitivity of the 
VSL calculations to various parameters that are 
employed to define the problem. Some of the 
parameters considered are the number of 
species in the chemistry model, surface catalytic 
activity, and freestream density. A variation in 
the values of these parameters gives an 
indication of the sensitivity of computed results 
to the uncertainty in their specified values. The 
sensitivity study has been done for altitudes of 
92.82, 88.45, and 59.60 km except for the 
chemistry model, which is studied by employing 
the freestream and surface conditions 
considered at the US-European Conference on 
high speed flowfields. 

Sensitivity to Number of Species in a 
Chemistry Model . Figure 4 shows the 
stagnation heating for a noncatalytic surface for 
two-temperature (2-T), 5-species (O, 02, N, N2, 
and NO), 7-species (O, O2, N, N2, NO, NO + , and 
e“), and 11-species (O, O2, N, N2, NO, 0 + , 02 + , 
N + , N2 + , NO+, and e") chemistry models. For 
the altitude (91.86 km) considered, surface and 


shock-slip are found to influence the surface 
heating more significantly than the number of 
species in the chemistry model. The effect on 
surface heating is negligible when the number of 
species is increased from 7 to 11. Therefore, a 
7-species chemistry model is considered 
adequate for the present study. 

Sensitivity to Surface Catalytic Activity. 

The surface heating calculations are made with 
the assumptions of noncatalytic, finite catalytic, 
and fully catalytic walls. The finite catalytic 
surface recombination values employed are 
those of Ref. 12. The VSL prediction using these 
values are compared with the STS-2 data 
(Ref. 19) in Fig. 5. For the 92.82 and 84.01 km 
altitude results shown in Figs. 6(a) and 6(b), the 
surface temperatures are less than about 
1100 K and the recombination probabilities for 
the finite catalytic wall in these figures are less 
than 3.25x10"® and 9.35X10" 4 for atomic oxygen 
and nitrogen, respectively. With these low 
values of the recombination probabilities, the 
calculated heat transfer rates are essentially the 
same as those for the noncatalytic wall case. 
However, at 59.60 km altitude (see Fig. 6(c)), the 
maximum recombination probability values for 
atomic oxygen and nitrogen increase to 
2.20x10'^ and 1.50x10'®, respectively, at the 
stagnation point where the temperature is about 
1460 K For a noncatalytic wall, a 34% 
reduction in stagnation VSL heating prediction 
is obtained for this case. However, for a fully 
catalytic wall this heating rate is increased by 
about 31%. For 92.82 km altitude (Fig. 6(a)) 
also, the VSL heating rates for the fully 
catalytic wall condition are larger, substantially 
so for the 84.01 km altitude case (Fig. 6(b)), as 
compared to the noncatalytic (or finite catalytic) 
wall values. It may be mentioned here that the 
VSL calculations at higher altitudes (Figs. 6(a) 
and 6(b)) employ a two-temperature (2-T) 
formulation with surface and shock-slip 
boundary condition and a one-temperature (1-T) 
model with no-slip boundary condition at lower 
altitudes (Fig. 6(c)) to reflect appropriate 
flowfield physics in the computations. 

Figures 6(a)-6(c) also contain the one- 
temperature (1-T) Navier-Stokes calculations of 
Yamamoto® for a noncatalytic wall with zero 
surface slip. Also shown in these figures are the 
OREX stagnation inferred flight data for 
comparison. Present VSL calculations for a 
finite-catalytic wall are in good agreement with 
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reported flight data as shown for altitudes of 
92.82 and 59.60 km (Fig. 6(a) and 6(c)). At 

84.01 km altitude (Fig. 6(b)), however, these 
calculations are lower by about 17% as 
compared to the data. Assuming the surface to 
be noncatalytic at lower altitudes will result in 
lower heating as compared to the inferred flight 
data as seen from Fig. 6(c). The Navier-Stokes 
calculations of Ref. 5 also give surface heating 
distributions, which appear at variance (see 
Figs. 6(a) and 6(c), for example) with the 
catalytic boundaiy condition employed in the 
present calculations. 

Sensitivity to Freestrea m Density . 
Actual atmospheric density during reentry can 
be a significant uncertainty for a given 
altitude/time, particularly at the higher 
altitudes. It may differ significantly from that 
obtained from a standard atmospheric model 
1*. The significance of this effect on surface 
heating is shown in the 'Comparison with Flight 
Data' section later. 

Assessment of Thermochemical 
Noneauilibrium and Slip Effects . An evaluation 
of thermochemical nonequilibrium in the 
flowfield can be made by analyzing the 
temperature profiles and surface distributions of 
pressure and heating rate for a non-catalytic 
wall. Figures 7(a) and 7(b) show stagnation 
profiles for the translational-rotational (T) and 
vibrational-electronic-electron (T ve ) 
temperature ratios for altitudes of 92.82 and 

84.01 km, respectively. At 92.82 km altitude 
(Fig. 7(a)), the two temperature ratios are quite 
different through the entire shock layer and 
accordingly, thermal nonequilibrium extends all 
the way to the surface. The degree of 
nonequilibrium in the two temperature profiles, 
however, is much reduced closer to the surface 
at 84.01 km altitude as shown in Fig. 7(b). 

Effects of thermal nonequilibrium on 
surface quantities may be evaluated by 
analyzing the pressure distributions given in 
Figs. 8(a) and 8(b). Unlike the continuum 
calculations, where assumptions are made with 
regard to the pressure tensor, a DSMC 
calculation accounts for nonisotropic effects 
which become significant under rarefied 
(thermal nonequilibrium) conditions. Therefore, 
normal force per unit area is shown in these 
figures from the DSMC calculations for 
comparison with the surface pressure (p w ) 
distribution obtained from the VSL 


computations. Under thermal equilibrium 
conditions, the isotropic pressure p w and normal 
force per unit area can be identified to be the 
same^. The surface distributions of these two 
quantities are different at the altitude of 
92.82 km (Fig. 8(a)) suggesting the influence of 
thermal nonequilibrium on surface pressure, 
p w , which is increased by about 20%. At 

84.01 km altitude (Fig. 8(b)), however, 
distributions of normal force/area and p w 
essentially have the same values, which imply 
conditions closer to thermal equilibrium. 

To assess the effect of thermal 
nonequilibrium on surface heating rate at 

84.01 km altitude, heating-rate distributions are 
computed (Kg. 9) using 1-T (for thermal 
equilibrium) and 2-T (for thermal 
nonequilibrium) models with the VSL method. 
Also, results are obtained with and without 
shock and surface-slip conditions. As shown in 
Kg. 9, the VSL results imply that the effects of 
both thermal nonequilibrium (even though 
present to some extent within the shock layer 
(Kg. 7b)) and slip are insignificant on surface 
calculated quantities for the 84.01 km altitude 
conditions and these effects may be neglected in 
calculations below this altitude. 

Comparison with Direct S imulation 
Monte Carlo (DSMC) and Navier- Stokes (NS) 
Calculations . Surface heating-rate distributions 
obtained using the VSL method are compared 
with those obtained from DSMC and NS 
calculations in Figs. 10(a) and 10(b) at altitudes 
of 92.82 and 84.01 km, respectively. Also shown 
are the flight inferred data for the stagnation- 
point heating. Except for the results shown over 
the conical flank in Kg. 10(a), there is generally 
a good agreement between the VSL and DSMC 
calculations for the two altitudes (Figs. 10(a) 
and 10(b)). This implies that a 2-T VSL 
formulation with slip boundary conditions 
models the low-density effects (thermal 
nonequilibrium and slip) quite well. The NS 
calculations of Ref. 5 are also in good agreement 
with the present results in the stagnation region 
at 92.82 km altitude (Kg. 10(a)) and along most 
parts of the OREX forebody at 84.01 km altitude 
(Fig. 10(b)). It is not obvious, however, why the 
NS results of Ref. 5 away from the stagnation 
region are higher in Fig. 10(a). Generally, the 
influence of a shoulder expansion on the 
upstream surface heating should become 
negligible at higher altitudes. The surface 
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heating becomes a function of only surface 
inclination at free molecular flow conditions. 
Present calculations agree quite well with the 
stagnation flight data at 92.82 km altitude. 
However, the present results are lower than the 
flight data by about 17% at 84.01 km altitude. 
Also, some slip effects are noticeable over the 
conical flank at this altitude 
(Fig. 10(b)). 

Additional comparisons of the present 
results with DSMC and NS calculations and 
flight data at the stagnation point are shown in 
the next section. 

Comparison with Flight Data 

This section provides a comparison 
between the computed and flight-inferred 
stagnation-point heating rates and electron 
number densities (measured from the 
electrostatic probe mounted near the shoulder of 
the vehicle). 

Stagnation-Point Heating . Figures 11(a) 
through 11(c) present VSL results (for altitudes 
of 105 to 48.4 km) as function of flight time from 
launch. The DSMC calculations for altitudes of 
105 to 79.9 km are shown in Fig. 11(c). Also 
shown in Figs. 11(a) and 11(c) are the NS 
calculations (called CFD-FEM results in Ref. 3) 
ofYamamoto®. 

The NS results of Refs. 3 and 5 were 
obtained with the thin-layer Navier-Stokes 
equations coupled with a finite-element 
modeling (FEM) solver that models the heat- 
transfer within the thermal protection material 
and provides the wall temperature boundary 
condition for the Navier-Stokes solver. The 
FEM solver accounted^ for both temperature 
and directional property dependence of the 
thermal protection materials considered. An 
outline of the coupling procedure between the 
NS and FEM solvers is given in Ref. 3. The NS 
results employed the following modeling 
features: no-slip and noncatalytic surface 
boundary conditions, and a one-temperature 
(1-T), 7-species nonequilibrium chemistry 
model. The unpublished results of Ref. 5 (used 
in the present work for comparison) were 
obtained recently and differ from those reported 
in Ref. 3 in that the agreement between the 
computed and flight-measured back surface 
temperature was improved for the reasons 
mentioned earlier. 

The inferred flight data shown in Fig. 11 
(and earlier) are taken from Ref. 3 and are 


referred to as ’inferred from flight’. The 
procedures used to obtain these data are 
described in Ref. 8. 

Since atmospheric data was not 
gathered during the OREX flight experiment, 
calculations of Refs. 3 and 5 are based on the 
atmospheric conditions given in Ref. 3, which is 
referred to as "OREX Atmosphere" in Figs. 11(a) 
and 11(b) and in Table 1. Also shown in these 
two figures are the VSL calculations based on 
the 'Jacchia Atmosphere’ model for altitudes 
above 90 km (see Table 1 for the two models). 

As can be seen from Fig. 11(a), the 2-T VSL 
predictions for a noncatalytic wall (with slip) 
based on the Jacchia atmosphere are in better 
agreement with inferred flight data between the 
90 and 105 km altitude range as compared to 
both the NS and VSL calculations utilizing the 
OREX atmospheric data. Further, the effect of 
slip on surface heating is dominant, whereas the 
thermal nonequilibrium (2-T) effects are 
secondary in this altitude range. Both slip and 
thermal nonequilibrium (2-T) effects are 
neglected in the NS 
calculations^* 5. 

For altitudes below 84.01 km, surface 
temperatures are higher and slip and thermal 
nonequilibrium effects become negligible. At 
higher temperatures, finite surface catalytic 
activity also begins to influence the heating 
rate. Therefore, a one-temperature (1-T), no- 
slip, and finite catalytic wall flow model is 
appropriate for the lower altitude calculations. 
The VSL calculations with these flow 
assumptions are in better agreement with the 
inferred flight data than the 1-T, no-slip, 
noncatalytic VSL and NS calculations as shown 
in Fig. 11(a). The dotted line in Fig. 11(b) is a 
fairing through the VSL predicted values with 
appropriate flow physics from high-to-low flight 
altitudes. As mentioned earlier, the low-density 
(i.e. slip and thermal nonequilibrium) effects 
become negligible for the OREX forebody 
surface results at an altitude of about 84 km. 

Figure 11(c) shows a comparison 
between VSL (with appropriate physics), NS, 
and DSMC predictions and with the inferred 
flight data. The DSMC results are shown for 
the high-altitude range of 105 to 79.9 km. 
Overall, the qualitative behavior of the DSMC 
results is similar to that of VSL calculations and 
there is a good agreement between the two and 
with the inferred flight data. The NS (1-T, no 
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slip, NCW) calculations of Yamamoto® are 
higher for altitudes greater than 88.45 km as 
compared to the VSL and DSMC calculations as 
well as the inferred flight data. As explained 
earlier, non-inclusion of thermal nonequilibrium 
and slip effects at higher altitudes and non- 
implementation of finite-catalytic wall boundary 
condition at lower altitudes with the NS 
computations^! 5 ma y ^ responsible for the 
differences with the VSL and DSMC predictions. 

As described in the Introduction section 
earlier, the OREX vehicle was enclosed in a 
protective fairing during launch, thereby 
negating the requirement for waterproofing the 
thermal protection materials. Since the Shuttle 
Orbiter, is an operational vehicle, waterproofing 
procedures are employed. Therefore, the OREX 
thermal protection system should not produce 
some of the outgasing products and consequent 
reduction in heating that is evident for the 
Orbiter^ during the initial portion of the heat 
pulse. The surface heat-transfer coefficient 
based on the initial entry heating data reported 
for OREX exhibits 1 ^ a monotonic increase with 
the increase in Knudsen number. This type of 
behavior is characteristic^, 21 0 f the 
hypersonic cold-wall, stagnation-point heating 
for a nonblowing surface. 

Electron Number Density . Flight data 
for the electron number density distribution in 
the boundary layer with strong entropy layer 
swallowing serve to delineate the applicability of 
f two-temperature (2-T) r and one-temperature 
(1-T)' flow models for nonequilibrium ionized gas 
flows. An accurate calculation of the electron 
density also provides knowledge of the onset and 
expiration of the communication black-out' 
phase during entry. 

Figures 12(a)- 12(c) show the electron 
number density history at positions where 
probes 1, 3, and 5 (see Fig. 1(c)) are located. 

Two sets of data shown in these figures are 
reduced^* 5 hy using 1-T and 2-T CFD methods 
of Ref. 6, which employs Navier-Stokes 
equations with a 7-species chemical model and 
Park’s 11 two-temperature model for the 2-T 
calculations. Also shown in these figures are 
predictions from these calculations reported in 
Ref. 2. Since the signal conditioner of the 
electrostatic probe is adjusted to measure the 
electron density in the 10 ^ - 10^0 particles/m^ 
range, there are no data measured for Probe 1 


between the flight times of 7401 and 7411 
seconds as shown in Fig. 12(a). It is not obvious, 
however, why the electron density first 
decreases at a flight time of about 7391 s and 
increases again at about 7421 s for Probe 1 
location. Similar decrease and increase in data 
is noticed for Probe 2 (not shown here) and 
Probe 3 (see Fig. 12(b)) also at different times in 
flight®. Figure 12(c) shows a continuous 
increase in electron density for Probe 5 with the 
increasing flight time (or decreasing altitude) as 
expected. 

Present VSL calculations are given at 
the three altitudes of 92.82, 88.45, and 84.01 km 
or the flight times of 7391, 7401, and 7411.5, 
respectively, in Figs. 12(a)- 12(c). These 
calculations are done with 1-T and 2-T flow 
models and slip boundary conditions. Present 

1- T VSL predictions are slightly higher than the 

2- T VSL results. Both of these VSL calculations 
are closer to the data and 1-T CFD calculations 
than the 2-T CFD results for all the three probe 
locations shown in Figs. 12(a)- 12(c). Similar to 
the surface heating results given earlier, the 
effect of thermal nonequilibrium on electron 
density also appears to be secondary. The 
observation of Ref. 2 that the flight data agree 
better with 1-T flow model is not necessarily the 
case for the present calculations. For flight 
times earlier than 7401 s (i.e. at higher 
altitudes), when the low-density effects (such as 
thermal nonequilibrium and slip) become 
important, a 2-T flow model with slip should be 
more realistic as evidenced by the good 
agreement between the 2-T VSL and DSMC 
calculations shown earlier. 

Finally, Figs. 13(a)- 13(c) show further 
detailed comparisons between the predicted 
values and data for the electron density profiles 
at the altitude of 84.01 km for the three probe 
locations. The VSL predictions in these figures 
have been obtained with the assumptions of 
thermal nonequilibrium (2-T) and equilibrium 
(1-T) with and without the slip boundary 
conditions. Except for the narrow region near 
the shock, the electron density profiles (similar 
to the surface heat-transfer rates) are predicted 
essentially by the thermal-equilibrium, no-slip 
flow model for the three probe locations. Except 
for Probe 1 (Fig. 13(a)), the electron density 
profiles pass through the flight-measured values 
for Probes 3 and 5 as shown in Figs. 13(b) and 
13(c). As noticed in Fig. 12(a) earlier, the reason 
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for the low flight values for Probe 1 as compared 
with the calculations is not obvious. Further, 
present VSL predictions are in good agreement 
with the NS results^ at this altitude for all the 
three probes. 

Summary and Conclusions 

This study presents a viscous shock- 
layer (VSL) analysis of the reentry flowfield 
around the forebody of the Japanese Orbital 
Reentry Experiment (OREX) vehicle. The 
OREX vehicle is a 50-deg half-angle spherically 
blunted cone with a nose radius of 1.35 m and a 
base diameter of 3.4 m. Obtained results span 
an altitude range of 105 to 48.4 km. In this 
altitude range, the flowfield character changes 
from thermal nonequilibrium to thermal 
equilibrium and the slip effects become 
insignificant at lower altitudes. The low density 
effects (such as thermal nonequilibrium and 
slip) at higher altitudes in the VSL method are 
accounted for through the two-temperature (2-T) 
formulation with slip boundary conditions. 

With the disappearance of these effects at lower 
altitudes, VSL results are computed from one- 
temperature (1-T) flow model with no-slip 
boundary conditions. Within the altitude range 
considered, the influence of low-density and the 
surface catalycity effects on flowfield structure 
is analyzed in detail. Present results are 
compared with the flight data and existing 
Navier-Stokes (NS) calculations for the OREX 
trajectory, and with the direct simulation Monte 
Carlo (DSMC) results at higher altitudes in this 
trajectory. 

Results from sensitivity studies indicate 
that a 7-species chemical model is adequate to 
analyze the OREX flowfield. However, the 
inclusion of low-density effects (i.e. both thermal 
nonequilibrium and slip) is important for 
altitudes greater than 84 km. This is 
demonstrated by the comparison of viscous 
shock-layer results obtained from one- 
temperature, no-slip and two- temperature, slip 
flow models. Surface catalytic activity also 
significantly influences the level of heating at 
altitudes lower than 84 km in the trajectory. At 
the peak heating altitude of 59.60 km, a 50% 
reduction in heating is obtained for a 
noncatalytic wall as compared to a fully 
catalytic wall. Reduction in heating is about 
23% for a finite-catalytic surface with 


recombination probabilities similar to those of 
the Shuttle Orbiter thermal protection coatings. 
However, at higher altitudes, the recombination 
probabilities are quite low and the calculated 
heating for altitudes greater than about 84 km 
is essentially the same as that for a noncatalytic 
surface. 

Since atmospheric data were not 
gathered in the OREX flight experiment, 
sensitivity of the atmospheric data used in 
earlier Navier-Stokes calculations is evaluated 
by employing the Jacchia atmospheric model for 
altitudes greater than 90 km. Generally, a 
better agreement between the VSL predictions 
and flight heating data is obtained using the 
Jacchia values. 

In general, there is a good agreement 
between the VSL predictions and the OREX 
measured or inferred data for electron number 
density and stagnation-point heat-transfer rate. 
Present results compare quite well with the 
direct simulation Monte Carlo (DSMC) 
predictions for altitudes higher than about 
84 km. The earlier Navier-Stokes calculations 
are higher at high altitudes due to the non- 
inclusion of low-density effects and lower at low 
altitudes due to the use of non-catalytic wall 
boundary conditions when compared with the 
flight-inferred heat-transfer rate data. 

The OREX heat-transfer rate data may 
be unique for high-altitude flight conditions in 
that they exhibit proper qualitative behavior 
with increasing altitudes. With further 
refinements in data extraction procedures and 
definitions, a valuable aerothermodynamic data 
base will become available. This data base 
should prove immensely helpful in validating 
the use of various theoretical models used for 
simulation and provide a test of flowfield codes 
for reacting gas flows. 
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Table 1. Flight Trajectory for OREX vehicle. 

(a) Free stream Conditions 


Flight Altitude, 

time, km 

s 

Velocity, 

m/s 

Too, 

K 

T w a 

K 



Mole Fractions 



Xo* 


X N; 

2 


Xo 

OREX 

Jacchia 

OREX 

Jacchia 

OREX 

Jacchia 

OREX Jacchia 


105.00 

7451.00 

218 

211.05 

332 


0.1528 


0.7815 

0.0 

0.0657 

7370.6 

101.10 

7454.65 

195 

196.89 

402 

0.2375 

0.1726 

0.7625 

0.7839 

0.0 

0.0435 

7381.0 

96.77 

7456.30 

192 

190.26 

485 

0.2375 

0.1884 

0.7625 

0.7863 

0.0 

0.0253 

7391.0 

92.82 

7454.10 

189 

188.30 

586 

0.2375 

0.2025 

0.7625 

0.7881 

0.0 

0.0094 

7401.0 

88.45 

7444.30 

187 


687 

0.2375 


0.7625 


0.0 


7411.5 

84.01 

7415.90 

189 


785 

0.2375 


0.7625 


0.0 


7421.5 

79.90 

7360.20 

199 


878 

0.2375 


0.7625 


0.0 


7431.5 

75.81 

7245.70 

207 


976 

0.2375 


0.7625 


0.0 


7441.5 

71.73 

7049.20 

215 


1091 

0.2375 


0.7625 


0.0 


7451.5 

67.66 

6720.30 

226 


1213 

0.2375 


0.7625 


0.0 


7461.5 

63.60 

6223.40 

237 


1344 

0.2375 


0.7625 


0.0 


7471.5 

59.60 

5561.60 

248 


1458 

0.2375 


0.7625 


0.0 


7481.5 

55.74 

4759.10 

259 


1531 

0.2375 


0.7625 


0.0 


7491.5 

51.99 

3873.40 

268 


1557 

0.2375 


0.7625 


0.0 


7501.5 

48.40 

3000.00 

271 


1388 

0.2375 


0.7625 


0.0 



a CFD inferred stagnation surface temperatures of Yamamoto (Ref. 5). Also, temperature distributions are specified 
from the same CFD computations. 


Table 1. Concluded. 

(b) Freestream Conditions 


Altitude, 

km 

Specific heat ratio, 

Too 

Speed of sound, 
m/s 

Mach No., 

Moo 

Mol. wt., 
kg/k-mole 

Density, 
kg/m ^ 

OREX 

Jacchia 

OREX 

Jacchia 

OREX 

Jacchia 

OREX 

Jacchia 

OREX Jacchia 

105.0 

1.40 

1.4108 

296.734 

298.202 

25.11 

24.986 

28.96 

27.838 

3.1400xl0* 7 2.3350xl0' 7 

101.1 

1.40 

1.4071 

280.989 

285.819 

26.53 

26.082 

28.96 

28.195 

5.7100xl0* 7 4.8341xl0' 7 

96.77 

1.40 

1.4041 

279.053 

279.323 

26.72 

26.694 

28.96 

28.466 

1.3810x10*® 9.3644x1 O' 7 

92.82 

1.40 

1.4015 

276.385 

276.492 

26.97 

26.960 

28.96 

28.701 

3.0090x10*® 1.9465x10*® 

88.45 

1.40 


275.002 


27.07 

28.96 



4.3060x10*® 

84.01 

1.40 


276.506 


26.82 

28.96 



1.0953xl0* 5 

79.90 

1.40 


274.430 


26.82 

28.96 



1.8455x10*® 

75.81 

1.40 


289.365 


25.04 

28.96 



3.6576xl0* 5 

71.73 

1.40 


295.069 


23.89 

28.96 



6.5184x10*® 

67.66 

1.40 


302.444 


22.22 

28.96 



1.2164x10*4 

63.60 

1.40 


309.776 


20.09 

28.96 



2.0594X10* 4 

59.60 

1.40 


316.900 


17.55 

28.96 



3.3131X10* 4 

55.74 

1.40 


323.528 


14.71 

28.96 



5.3150x10*^ 

51.99 

1.40 


328.254 


11.80 

28.96 



8.2445x10^ 

48.40 

1.40 


331.126 


9.06 

28.96 



1.2677x10*3 


11 





1508 


14 


b 


NoncataJytic surface 



Fig 4. Effect of chemistry model and slip on stagnation 
heating with thermal nonequilibrium (2-T). 


Altitude = 92.35 km 
V M = 7.50 km/s 


MW/m 2 



Fig. 5 Comparison of VSL surface heating predictions 
with STS-2 flight data along the windward center-line plane. 


q, W/cm 
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o Non catalytic. NS, 1 -T, No Slip (Yamamoto) 
• Inferred flight data 
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(a) Alt = 92.82 km (Kn^ = 0.0086) 


2.5 


Fig. 6 Effect of surface catalytic activity on 
surface heating rate. 


Non calalytic, VSL, 2-T, Slip 



Fig. 6 Continued. 



Fig. 6 Concluded. 
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(a) Alt = 92.82 km (Kn„ = 0.0086) 


Fig. 7 Stagnation temperature profiles (VSL, 2-T, slip). 



(b) Alt = 84.01 km (Kn^ = 0.0016) 


Fig. 7 Conluded. 



s, m 

(a) Alt = 92.82 km (Kn^ = 0.0086) 
Fig.8 Surface pressure distribution. 



s, m 

(b) Alt = 84.01 km (Kn M = 0.0016) 
Fig. 8 Conluded. 



s, m 

Fig. 9 Predicted surface heating-rate distribution from 
1-T and 2-T VSL calculations with and 
without slip (Alt = 84.01). 


Nan Catalytic wall 


DSMC 



s, m 


(a) Alt = 92.82 km (Kn^ = 0.0086) 

Fig. 10 Comparison of surface heating-rate distribution. 



(b) Alt = 84.01 km (Kn M = 0.0016) 
Fig. 10 Concluded. 
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(a) VSL and NS calculations and flight data 

Fig. 1 1 Comparison of predicted stagnation-point heat 
transfer rates with inferred flight data. 
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Inferred from flight 

VSL calculations 

(with appropnate physics) 

VSL Predictions: 

OREX Atmosphere 
o 1-T, no slip, NCW 
$ 1-T, no slip, FiCW 

JaccNa Atmosphere 
* 2-T. slip, NCW 


7300 7350 7400 7450 7500 7550 
Time, s 

(b) Flight data and VSL results 
Fig. 1 1 Continued. 
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(c) VSL, NS, and DSMC calculations and flight data 
Fig. 11 Concluded. 
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LoQlO (electron number density (/m 3 )] Log^o (electron number density (/m 3 )J LoQtO (electron number density (/m 3 )] 


A Right data anatyzed by 1 -temp. CFD 



100 95 90 85 80 75 70 

Altitude (km) 

(a) Probe 1 (s = 1.8327 m, n = 1.1665 x 10‘ 2 m) 
Fig. 12. Electron number density history. 

* Flight data analyzed by 1 -temp. CFD 



7380 7400 7420 7440 

Elapsed time after Irft-off (sec) 

l - J L -J 1 I I I I I L -J- 1 

100 95 90 85 80 75 70 

Attitude (km) 

(b) Probe 3 (s = 1.8731 m, n = 3.50 x 10’^ m) 

Fig. 12. Continued. 


A Right data analyzed by 1-temp. CFD 
0 Right data analyzed by 2-temp. CFD 



Altitude (km) 

(c) Probe 5 (s = 1.9135 m,n = 5.8335 x 10"^ m) 
Fig. 12. Concluded. 



(a) Probe 1 (s = 1.8327 m, n = 1.1665 x 10 ^ m) 
Fig. 13 Comparison of electron number-density profiles 



(b) Probe 3 (s = 1.8731 m, n = 3.50 x 10"^ m) 
Fig. 13 Continued. 



n, m 

(c) Probe 5 (s = 1 .9135 m, n = 5.8335 x 10"^ m) 
Fig. 13 Concluded. 
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